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Abstract 

m 

Neutrinos interacting with the quintessence field can trigger the accelerated expansion 
of the Universe. In such models with a growing neutrino mass the homogeneous cosmolog- 

^vq ', ical solution is often unstable to perturbations. We present static, spherically symmetric 

solutions of the Einstein equations in the same models. They describe astophysical objects 
composed of neutrinos, held together by gravity and the attractive force mediated by the 
quintessence field. We discuss their characteristics as a function of the present neutrino 
mass. We suggest that these objects are the likely outcome of the growth of cosmological 
perturbations. 
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The mechanism responsible for the onset of the accelerating phase in quintessence 
cosmology remains undetermined. Explaining the emergence of an accelerating phase in 
recent cosmological times constitutes one of the most difficult challenges of quintessence 
models - the coincidence problem. A possible trigger for the acceleration has been pro- 
posed recently [1, 2], arising through the interaction of the quintessence field with a matter 
component whose mass grows with time. This matter component may be identified with 
neutrinos [1, 2, 3]. In the proposed scenario the neutrinos remain essentially massless until 
recent times. When their mass eventually grows close to its present value, their interaction 
with the quintessence field (the cosmon) almost stops its evolution. The potential energy 
of the cosmon becomes the dominant contribution to the energy density of the Universe. 
Cosmological acceleration ensues. 

For the coupled neutrino-cosmon fluid the squared sound speed c 2 may become nega- 
tive - a signal of instability [3]. Indeed, the sign of c 2 oscillates in the accelerating phase 
for one of the proposed models [2]. A natural interpretation of this instability is that the 
Universe becomes inhomogeneous with the neutrinos forming denser structures. Within 
the linear approximation the neutrino fluctuations can be followed in these models until a 
redshift around one, when the neutrino overdensities become nonlinear [4]. One suspects 
that some form of subsequent collapse of these fluctuations will result into bound neutrino 
lumps. In this letter we present static, spherically symmetric solutions of the Einstein 
equations that describe such structures and study their characteristics. Astrophysical 
objects composed of neutrinos have also been studied in [5, 6]. 

We assume that the energy density of the Universe involves a gas of weakly interacting 
particles (neutrinos). The mass m of the particles depends on the value of a slowly varying 
cosmon field (f> [7]. For the field equation 
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we approximate the neutrino energy- momentum tensor as T^ v = diag(— p, p, p,p). The 
cosmology of [1, 2] also assumes the presence of another gas of particles (dark matter) 
whose mass is independent of <j>. 

We consider stationary, spherically symmetric configurations, with metric 

ds 2 = -B{r)dt 2 + r 2 (d6 2 + sin 2 9 dip 2 ) + A(r)dr 2 . (2) 

For the neutrinos we assume a Fermi-Dirac distribution, with locally varying density - the 
Thomas-Fermi approximation. The local chemical potential satisfies fi(r) = n®/ \/B(r) 
[6, 8]. Stable configurations are prevented from collapsing by the pressure generated 
through the exclusion principle. We concentrate on vanishing temperature of the neutrino 
gas. We do not expect qualitative changes of our solution for a non-zero temperature. For 
simplicity we consider one neutrino species, with the generalization to degenerate neutrino 
masses being straightforward. 

We parametrize the particle mass by a dimensionless function in, defined according 
to m((f>) = am [((/) — cj))/M], with a an arbitrary energy scale and M = (16-7rG) -1 ' 2 ~ 
1.72 x 10 18 GeV. Here is a fixed reference value, close to the present value of the 
quintessence field. Hence, a is of the order of the present neutrino mass, in the eV 
range or somewhat below. For concreteness, we consider a cosmon potential of the form 



U{4>) = Ca 4 exp [— a(<p — 4>)/M] with a = 0(1). However, the effect of the potential 
on our solutions is negligible. For this reason, the predicted astrophysical objects are 
largely independent of the form of the potential, and depend mainly on the interaction 
between dark energy and neutrinos. The present cosmological value of (f> is given by the 
requirement that U(<j>) constitute about 3/4 of the critical energy density U(<j>) — 10 -11 
(eV) 4 . The cosmological value of <f) is taken as the asymptotic value (p as of our local 
solutions for large r, obeying (<p as — $)/M = <j) as ~ (l/o) [25.3 + InC + 4 In (cr/eV)]. 

The equations of motion become more transparent if we define the dimensionless vari- 
ables <fi = (<fi — 4>)/M and f = a 2 r/M. All other dimensionful quantities are multiplied 
with appropriate powers of a, in order to form dimensionless quantities denoted as tilded. 

We use B = B/jxq = B<j 2 //j,q and jtff) = 1/J B{r). We define the radius R of the compact 
object by the value of f at which the fermionic density becomes negligible. The physical 
radius is 

i?/Mpc ~ 1.1 x 1(T 2 (a/eV)~ 2 R. 

The mass of the object is given by its Schwarzschild radius R s . For f — ► oo we have 
B = 1/A = 1 — R s /f. In units of the solar mass, the mass of the neutrino lump is 



-2 f, 



Mtot/MQ ~ 1.2 x 10 1 ' (a/eV)-'R s . 



Another important characteristic is the total neutrino number, which we assume to be 
conserved. It is 

N si 5.1 x 10 81 (a/eVr 3 iV, 



with 



The field equations read [8] 
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where a prime denotes a derivative with respect to f. We also have 
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for p, > fn, and n = p = p = for jX < fh. Finally, U{4>) = Cexp ( —a<p). 
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Figure 1: Radial dependence of (J), p, A and B for a neutrino lump. 



We need four initial conditions for the system of equations (3). Two of them are 
imposed by the regularity of the solution at f = 0: 4>'(0) = 0, A(0) = 1. The value 
of -B(O) is the only free integration constant. Since AB(r — ► oo) = 1 one has AB(r — ► 
oo) = (fi /a)~ 2 . As a result, the choice of -B(O) determines the chemical potential and, 
therefore, the total number of neutrinos in the lump. Finally, (j)(0) must be chosen so that 
4>{f — ► oo) reproduces correctly the present value <j) as of the cosmological solution. (We 
assume that the time scale of the cosmological solution is very large and neglect the time 
dependence of (p(f — ► oo).) 

We consider two types of models, distinguished by the dependence of the particle mass 
on the field: 

Model I assumes m(</>) = —1/0 [2], with the field </> taking negative values. 
Model II assumes m(0) = exp (—&</>) [1], with b < 0. (Notice that a = a/y/2, b = (3/-\/2 
in comparison to [1], where a different convention for M is used.) 

In both cases we are interested in values of the field near <fi = 0. For model II we can 
choose 4> such that <fi as = 0, implying that a = m u (to) equals the present neutrino mass. 
One infers for the quintessence potential InC = —25.3 — 41n(o"/eV). The parameter a 
is fixed by requiring that during the early stages of the cosmological evolution the dark 
energy be subleading and track the radiation or the dark matter. During the radiation 
and matter dominated epochs, the dark energy follows a "tracker" solution with a con- 



stant density parameter Qh,eariy = ^/(^o 2 ), where n = 3 (4) for matter (radiation) [1]. 
Observations require a to be large, typically a ^ 7 [9]. We use a = 7 in the following. 
The future of our Universe is described by a different attractor, for which the dark energy 
dominates. Our present era coincides with the transition between the two cosmic attrac- 
tors. The influence of the neutrinos on the evolution of the cosmon field is determined 
by the second term in the r.h.s. of eq. (1). Demanding that today this term be equal 
to the first term, that arises from the potential, fixes the present neutrino fraction to 
the value £l u (to) = —(b/a)Q.h{to) [1]. For a realistic cosmology with present dark energy 
fraction f^(to) — 3/4 one has to adjust b to the neutrino mass. For one dominant neu- 
trino species we have b = —a(36eV/m u (to)). For model I we need to know how close 
(j) as is to zero, with a = —(^as'^uito). As compared to model II, we have now an effective 
^-dependent b(cj>) = —l/4>, which results in the condition cp as = — (l/o)m w (to)/(36eV) or 
a = (l/a)m 2 (i )/(36eV). 

In fig. 1 we present a typical solution describing a static astrophysical object in model 
I. The chemical potential has the value jXq ~ 2.9. The scalar field becomes more negative 
near the center of the solution, so that the neutrinos become lighter there. The asymptotic 
value is (p as = —0.02, which corresponds to m u (to) ~ 5 eV. The pressure and density of 
the fermionic gas vanish for f > R ~ 0.91. The mass of the object can be deduced from 
the asymptotic form of A or B for f — » oo. We find R s ~ 0.12. The total fermionic 
number is iV ~ 0.88. The form of the solutions in model II is similar to the one depicted 
in fig. 1. 

The variation of the chemical potential results in a whole class of solutions, depicted 
by the solid line in fig. 2. We display the dimensionless Schwarzschild radius R s as a 
function of the dimensionless radius of the object R. There is a maximal value for the 
mass, denoted by the end of the thick line. The continuation of the curve has the form of 
a spiral and is depicted by a thinner line. This branch is unstable to perturbations that 
can lead to gravitational collapse [10]. In order to demonstrate this fact, we plot in the 
same figure R s as a function of N/6 (dotted line). This curve has two branches. The one 
depicted by a thinner line corresponds to the thinner line of the curve R S (R). There are 
two possible values of R s that correspond to the same value of the total neutrino number 
N. The value on the thinner line has a larger value of R s and results in a larger mass. 
The corresponding configuration is unstable towards one with the same iV located on the 
thicker line. The characteristics of the solutions depend only very mildly on the value of 
4>as, as demonstrated by the comparison of the solid and dashed curves. All the values 
of m u (to) in the range [0,5] eV correspond to (p as in the range [-0.02,0]. The respective 
R S (R) curves lie between the solid and dashed curves of fig. 2. 

A striking feature is the existence of neutrino lumps with arbitrarily small mass. They 
correspond to the lower left corner of the figure, where both R and R s vanish. For such 
objects the contribution from gravity is negligible and their existence is a consequence of 
the attractive force mediated by the scalar field. Such configurations are not generic, but 
depend crucially on the assumed form of rh((p). A completely different form of solutions 
appears in model II. In fig. 2 we also depict the gravitational potential <3?(f) = —R s /{2f) 
at a distance f = R equal to the radius of the astrophysical object. 

The function R S (R) in model II displays a different behaviour. In fig. 3 we plot this 
function for four different values of 6, namely b = —500, —50, —4, — 1. Realistic neutrino 
masses correspond to large, negative b. For b = — 1 (solid line) there is a maximal 
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Figure 2: Mass vs. size for neutrino lumps in model I. 



value for the mass of the astrophysical objects and a branch of unstable solutions. The 
maximal value of R s is comparable for model I and model II with b = — 1, even though the 
corresponding radius is larger by an order of magnitude in the second case. For b = —4 
(dashed line) the maximal value of R s and the corresponding R increase by roughly two 
orders of magnitude. 

The crucial qualitative difference with model I concerns the form of the solutions 
with low values of R s . In model I for R s — > we have R — > 0, while in model II we 
have R — > oo. The attractive interaction mediated by the scalar field in model II is not 
sufficiently strong to lead to bound objects with a small fermion number. Gravity must 
play a role for compact objects to exist. As \b\ increases the dependence ofmoiKji becomes 
more pronounced. The effective neutrino mass in the interior of a compact object can 
become smaller without a large variation of c/> (and a significant energy cost through the 
field derivative term). This has two significant effects: a) Objects with smaller iV and R 
can exist. As a result the bending of the curve R S (R) for low R s takes place for smaller 
R. b) The configurations that are gravitationally unstable (indicated by the spiral in the 
upper part of the curve) are shifted toward larger values of R s . The reason is that the 
neutrinos are essentially massless in the interior of of such configurations, carrying only 
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Figure 3: Same for model II. 



kinetic energy. This makes the gravitional collapse difficult. 

The curves R S (R) in model II with b = —500 and —50 are also depicted in fig. 3. 
We have not managed to determine numerically a maximal value of R s , as objects with 
huge values of R s , R (larger by more than twenty orders of magnitude than the ones 
depicted) are possible. For comparison we note than in model I we have a maximal value 
(R s )max = 0.80 with a corresponding radius (R) ma x = 2.0. In fig. 3 we observe minimal 
values of the radius, (R) m i n = 0.54 for b = —50 and {R) m in = 0.054 for b = —500. 
The corresponding values of the Schwarzschild radius are (R s )min = 1-1 x 10~ 5 and 
(R s )min = 1-1 x 10 -8 , respectively It is apparent that for small R s we have the scaling 
behaviour R ~ b~ , R s ~ 6~ 3 . This can be understood by noticing that in the limit 
A', B' — ► 0, A — ► 1, and for negligible dU/dcj), the factors of b in eq. (3) can be eliminated 
through the redefinitions b(j) — > eft, br — v f. In fig. 3 we also depict the surface gravitational 
potential &{R) = —R S /{2R) as a function of R for the cases b = —500 and —50. 

In fig. 4 we display the size R of the astrophysical objects as a function of the present 
neutrino mass m u = m v (to). Restoring physical units requires the scale a. We use for 
model I a = 7 or (a/eV) 1 ' 2 ~ 0.063(m i/ /eV). The function R S (R) has a very mild 
dependence on <p as for < <f) as < 0.02 (see fig. 2). For given R the variation of a (or 
equivalently m u ) produces a class of astrophysical objects of variable physical size. They 
all generate the same surface gravitational potential <3? = —R S /(2R). In fig. 4 we depict 
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Figure 4: Size of neutrino lumps as a function of the neutrino mass. 



three such classes. The first two contain objects with strong gravitational potentials, while 
the last one contains objects that generate weaker fields. Solutions with R — > produce 
curves parallel to those in fig. 4, but located closer to the lower left corner. In the same 
figure we also depict two solutions of model II. In this model the neutrino mass is uniquely 
determined by the value of b. The two points in fig. 4 correspond to the minimal values 
of R for b = —50 and —500. These are R = 0.54 and R = 0.054, respectively. 

Recently, a first investigation of the coupled fluctuations of dark matter, neutrinos, 
baryons and the cosmon field has been performed for the models within the linear ap- 
proximation [4]. For a specific model with a present average neutrino mass of 2.1 eV, the 
neutrino fluctuations grow nonlinear at a redshift around one. The typical size of these 
fluctuations is large, in the range of superclusters and beyond. A further investigation 
of the fate of these neutrino lumps will have to follow their collapse due to the scalar- 
mediated attractive interaction and gravity. This should generate the distribution of the 
integration constants of the present solution, like the characteristic mass and size of the 
lumps. 

Our study demonstrates that the presence of instabilities in quintessence cosmologies 
with a variable neutrino mass may have interesting astrophysical consequences. After 
a sufficiently long time, these instabilities may lead to the formation of stable bound 
neutrino lumps. Their radius and mass within the family of allowed solutions (for given 
vfiy) depend on the details of the dynamical formation mechanism. Since in the models of 
[1, 2] the neutrinos remain free streaming until a rather recent cosmological epoch (say, 



z = 5), one may expect a large typical size of the neutrino lumps (more than 100 Mpc). At 
the present stage of the investigations it is not clear if such lumps have already decoupled 
from the cosmological expansion - for this, the perturbations have to grow nonlinear - 
or if this will happen only in the future. In the extreme case of an early formation of a 
population of lumps with subgalactic size, they could even play the role of dark matter. 
The detection of lumps could proceed directly through their gravitational potential, or 
indirectly through their attraction for baryons. Quintessence cosmologies may provide 
surprises for structures on very large scales. 
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